The final stages of accretion onto non-Kerr compact objects 



Cosimo Bambi 1,2 'E] and Enrico Barausse 3 '[l] 

1 Institute for the Physics and Mathematics of the Universe, 
The University of Tokyo, Kashiwa, Chiba 277-8583, Japan 
2 Arnold Sommerfeld Center for Theoretical Physics, 
Ludwig- Maximilians- Universitdt Miinchen, 80333 Munich, Germany 
3 Maryland Center for Fundamental Physics & Joint Space-Science Institute, 
Department of Physics, University of Maryland, College Park, MD 20742, USA 

(Dated: October 11, 2011) 

The 5 — 20 Mq dark objects in X-ray binary systems and the 10 5 — 1O 9 M0 dark objects in galactic 
nuclei are currently thought to be the Kerr black holes predicted by General Relativity. However, 
direct observational evidence for this identification is still elusive, and the only viable approach to 
confirm the Kerr black hole hypothesis is to explore and rule out any other possibility. Here we 
investigate the final stages of the accretion process onto generic compact objects. While for Kerr 
black holes and for more oblate bodies the accreting gas reaches the innermost stable circular orbit 
(ISCO) and plunges into the compact object, we find that for more prolate bodies several scenarios 
are possible, depending on the spacetime geometry. In particular, we find examples in which the 
gas reaches the ISCO, but then gets trapped between the ISCO and the compact object. In this 
situation, accretion onto the compact object is possible only if the gas loses additional angular 
momentum, forming torus-like structures inside the ISCO. 

PACS numbers: 04.50.Kd, 97.10.Gz, 97.60.Lf 



I. INTRODUCTION 

There is by now robust observational evidence sup- 
porting the existence of 5 — 2OM dark objects in X-ray 
binary systems and that of 10 5 — 1O 9 M dark objects in 
galactic nuclei [1 . These objects are thought to be the 
black holes (BHs) predicted by General Relativity (GR), 
since their characteristics cannot be explained otherwise 
without introducing new physics. In 4-dimensional GR, 
an (uncharged) stationary BH is described by the Kerr 
solution and is uniquely characterized by two parame- 
ters, the mass M and the spin angular momentum J [2]. 
A fundamental limit for a Kerr BH is the bound \a\ < 1, 
where a = cJ/(GM 2 ) is the spin parameter, which en- 
sures that the Kerr metric describes a BH and not a 
naked singularity. 

However, the evidence that the geometry around these 
BH candidates is described by the Kerr metric is still cir- 
cumstantial, and the Kerr BH hypothesis mainly relies 
on the assumption that GR is the correct theory of grav- 
ity. Because GR has been tested only in the weak-field 
regime [3], several authors have suggested possible ways 
to further test the Kerr BH hypothesis using future, and 
in some cases even present, data. 

More specifically, the detection of extreme mass ratio 
inspirals (EMRIs, i.e. systems consisting of a stellar- 
mass compact object orbiting a supermassive BH candi- 
date) with future space-based gravitational-wave detec- 
tors will allow one to map the spacetime geometry with 
exquisite accuracy. This is because missions like LISA 
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or a similar European mission will be able to follow the 
stellar-mass compact object for millions of orbits around 
the central supermassive BH candidate, and therefore de- 
viations from the Kerr geometry will lead to a phase dif- 
ference in the gravitational waveforms that grows with 
the number of observed cycles. This technique is very 
promising and has been studied in detail by many au- 
thors [4H6] . Likewise, future gravitational- wave detec- 
tors will be able to detect the quasi-normal modes of BH 
candidates, and because for a Kerr BH the frequencies of 
these modes depend only on the spacetime's mass M and 
spin J, any departure from this pattern will allow one to 
measure deviations away from the Kerr geometry [7 . 

The geometry around BH candidates can be con- 
strained also with present or future electromagnetic data. 
For instance, radio timing observations allow constraints 
to be put on the quadrupole moment of the compact com- 
panions of radio pulsars [8 , while astrometric monitoring 
of stars orbiting at milliparsec distances from Sgr A* may 
be used to constrain models for the supermassive BH can- 
didate at the center of the Galaxy [9] . Constraints on the 
nature of BH candidates can also be obtained by extend- 
ing the methods currently used to estimate the spins of 
these objects, such as X-ray continuum [10] and iron-Ka 
measurements [11] , observations of quasi-periodic oscilla- 
tions [12 , and measurements of the cosmic X-ray back- 
ground [13] . These methods can in principle be applied 
even with present data, provided that the systematic er- 
rors are properly understood. Future observations of the 
shadow of nearby supermassive BH candidates are an- 
other exciting possibility to test the BH paradigm [14]. 

BH candidates are often surrounded by an accretion 
disk. For the stellar-mass objects in X-ray binary sys- 
tems, the disk originates from the material stripped from 
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the companion, while the gas accreting onto the super- 
massive objects in galactic nuclei comes from the inter- 
stellar medium. When the mass accretion rate is moder- 
ate and the accreting gas has significant angular momen- 
tum, the disk is geometrically thin and optically thick. 
The standard framework to describe these disks is the 
Novikov-Thorne model [15], where the disk lies on the 
equatorial plane of the system and the gas moves on 
nearly geodesic circular orbits. If the central object is 
a BH, the inner edge of the disk is assumed to be at 
the radius of the innermost stable circular orbit (ISCO): 
circular orbits inside the ISCO are radially unstable, so 
the gas reaches the ISCO moving on quasi- circular or- 
bits, and then quickly plunges into the BH. A crucial in- 
gredient of the Novikov-Thorne model is the assumption 
that the gas does not emit additional radiation as soon 
as it enters the plunging region. Numerical simulations 
show that deviations from this picture are small relative 
to the effect of the uncertainties in other parameters of 
the system [16] and can therefore be neglected. (Note 
however that the authors of Ref. [17 reach a different 
conclusion, because using general relativist ic magneto- 
hydrodynamics simulations they find that there can be 
significant magnetic stress inside the ISCO.) 

If the central object is not a Kerr BH, the final stages 
of the accretion process may be different. In this work 
we focus on the so-called Manko-Novikov (MN) space- 
times [18]. These are stationary, axisymmetric, and 
asymptotically flat exact solutions of the vacuum Ein- 
stein equations with arbitrary mass-mult ipole moments, 
and they can therefore describe the geometry outside a 
generic compact object, be it a Kerr black hole or some 
other exotic object within GR. In particular, we con- 
sider a subclass of MN spacetimes characterized by three 
parameters (mass M, spin angular momentum J and 
mass quadrupole moment Q). In a Kerr spacetime, the 
quadrupole moment Qxerr and all the higher order mul- 
tipole moments are known to be functions of the mass 
and spin. Thus, if Q = Qxerr(^ J) our MN spacetimes 
exactly reduce to a Kerr BH. This property is very con- 
venient because it makes these spacetimes an ideal tool 
to set-up a null experiment to test the BH paradigm: any 
experiment pointing at a value of Q significantly differ- 
ent from QKerr(^ J) would imply that the object under 
consideration is either a compact object different from a 
BH within GR, or a BH or a compact object in a gravity 
theory different from GR. In this sense, MN spacetimes 
can be used not only to test the BH hypothesis but also 
to test the gravity theory itself. 



The accretion process in our MN spacetimes can be 
qualitatively different than in a Kerr spacetime. In par- 
ticular, we find that when the accreting gas reaches the 
ISCO (i.e. the inner edge of the Novikov-Thorne disk 
model) there are four qualitatively different possibilities: 

1. The ISCO is radially unstable, and the gas plunges 
into the compact object remaining roughly on the 
equatorial plane. This is the same scenario as in 
the Kerr case. 

2. The ISCO is radially unstable and the gas plunges, 
but does not reach the compact object. Instead, it 
gets trapped between the object and the ISCO and 
forms a thick disk. 

3. The ISCO is vertically unstable, and the gas 
plunges into the compact object outside the equa- 
torial plane. 

4. The ISCO is vertically unstable and the gas 
plunges, but does not reach the compact object. 
Instead, it gets trapped between the object and the 
ISCO and forms two thick disks, above and below 
the equatorial plane. 

This paper is organized as follows. In Sec. [II] we review 
a subclass of the MN solutions, characterized by three 
free parameters (M, J, Q), that we use to describe the 
spacetime around generic compact objects. In Sec. |III[ 
we study plunging orbits in these MN spacetimes and we 
show that the accreting gas may not reach the surface of 
the compact object. This may lead to the formation of 
a thick disk inside the ISCO. In Sec. IIVI we discuss the 
possible astrophysical consequences of these thick disks. 
Finally, in Sec. [V] we draw our conclusions, while in Ap- 
pendix [A] we review the theory of thick non-gravitating 
disks in axisymmetric stationary spacetimes, which we 
use to describe the thick disks inside the ISCO. In Ap- 
pendix [B] we review the thermal bremsstrahlung emission 
rate. 

Throughout the paper we use units in which G = c = 1 , 
unless stated otherwise. 



II. MANKO-NOVIKOV SPACETIMES 

The MN metric [18] is a stationary, axisymmetric, 
and asymptotically flat exact solution of the Einstein 
vacuum equations with arbitrary mass-multipole mo- 
ments. In quasi-cylindrical coordinates (p, z) and in pro- 
late spheroidal coordinates (x,?/), the line element is, re- 
spectively, 
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Here R = \J x 2 + y 2 — 1 and P n are the Legendre poly- 
nomials with argument xy/R, 
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while the relation between prolate spheroidal and quasi- 
cylindrical coordinates is given by 
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The MN solution has an infinite number of free param- 
eters: fc, which regulates the mass of the spacetime; a, 
which regulates the spin; and a n (n = 1, ...,+oo) which 
regulates the mass-multipole moments, starting from the 
dipole ai, to the quadrupole «2, etc. For a ^ and 
a n = 0, the MN solution reduces to the Kerr metric. For 
a = a n = 0, it reduces to the Schwarzschild solution. For 
a = and a n ^ 0, one obtains the static Weyl metric. 

The no-hair theorem [19] states that the only asymp- 
totically flat, vacuum and stationary solution of the Ein- 
stein equations that is non-singular on and outside an 
event horizon and that presents no closed timelike curves 
outside it is given by the Kerr metric. Therefore, the MN 



spacetime must either have no event horizon, or present 
naked singularities or closed timelike curves outside it. 
In fact, the surface x = 1 (p = 0, \z\ < fc), which is the 
event horizon in the Kerr case a n = 0, in general is only 
a partial horizon, because it presents a naked curvature 
singularity on the equatorial plane (i.e. at x = 1, y = 0, 
corresponding to p = z = 0) [T5] . Also, there are closed 
timelike curves outside it. However, these pathological 
features appear at small radii and here the basic idea is 
that naked singularities and closed timelike curves do not 
exist in reality because they are either inside some sort of 
exotic compact object, whose exterior gravitational field 
is described by the MN metric, or because GR breaks 
down close to them; see e.g. Ref. [20 for some specific 
mechanisms that can do the job. 

Without loss of generality, we can put ql\ = to bring 
the massive object to the origin of the coordinate system. 
In what follows, we restrict our attention to the subclass 
of MN spacetimes with a n = for n ^ 2. We then have 
three free parameters (fc, a, and 0L2) related to the mass 
M, the dimensionless spin parameter a = J/M 2 , and 
the dimensionless anomalous quadrupole moment q = 
— (Q — QKerr)/^ 3 , by the relations 



vT 



k = M- 



1 



a 



a 2 



1 + a 2 
1 fc 3 ' 



(2.13) 
(2.14) 
(2.15) 



4 



Note that q measures the deviation from the quadrupole 
moment of a Kerr BH. In particular, since Qxerr = 
— a 2 M 3 , the solution is oblate for q > —a 2 and prolate 
for q < —a 2 . When q = 0, the solution reduces to the 
Kerr metric, but when q ^ also the higher-order mass- 
mult ipole moments have a different value than in Kerr. 

Because the MN metric is stationary and axisymmet- 
ric, geodesic orbits have two constants of motion, the spe- 
cific energy E = —u t and the z-component of the specific 
angular momentum L = u^. The t- and ^-components 
of the 4-velocity of a test-particle are therefore 
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From the normalization of the 4-velocity, g^ v u^u v — — 1, 
we can write 
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where the effective potential V e s is defined by 
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Circular orbits in the equatorial plane must have p = 
z = 0, which implies V e s = 0, and p = z = 0, 
which implies d p V e ft = and d z V e f? = 0. This means 
that circular equatorial orbits are located at simultane- 
ous zeros and extrema of the effective potential, where 
dpVes = d z V G f[ = Kff = 0. Because d z V e f? = is satis- 
fied identically for z = (simply because of the reflection 
symmetry of the MN metric with respect to the equato- 
rial plane), from these conditions one can obtain E and 
L as a function of the radius r of the circular equatorial 
orbit. These orbits are stable under small perturbations 
in the radial direction if d 2 V e R < 0, and in the vertical 
direction if d 2 V e ^ < 0. 

As far as the position of the ISCO is concerned, for 
any given value of a there are two critical values q\ and 
q2 for the anomalous quadrupole moment (both negative, 
Q2 < Qi < 0, but whose exact values depend on the spin 
a). For q > qi, circular orbits on the equatorial plane 
are always vertically stable and the ISCO radius is de- 
termined by the onset of the orbital instability in the 
radial direction (i.e. the ISCO is the marginally stable 
circular orbit in the radial direction) . This is the picture 
that one is familiar with in the Kerr case, which is indeed 
included in this range of g's (since for Kerr q = 0). 

For q 2 < q < qi, there are two circular orbits r = 
ri and r — r 2 (with r\ > r 2 ) that are vertically stable 
but only marginally stable in the radial direction, and 
one circular orbit r = rs (with rs < r 2 ) that is radially 
stable but only marginally stable in the vertical direction. 
Stable circular orbits therefore exist for r > r\ and for 
rs < r < r2, while orbits with r 2 < r < r\ are radially 
unstable (although vertically stable) and orbits with r < 




FIG. 1. The character of the effective ISCO for MN space- 
times with spin parameter a and anomalous quadrupole mo- 
ment q. Compact objects to the right (left) of the solid red 
line have ISCOs determined by the instability in the radial 
(vertical) direction. Compact objects to the right (left) of the 
dotted blue line are oblate (prolate) bodies, while they are 
more oblate (prolate) than a BH if q > (q < 0). 



rs are vertically unstable (although radially stable). As 
far as an accretion disk is concerned, however, what is 
relevant is the radius r = r\ of the outer marginally stable 
circular orbit, because that is the radius at which the gas 
starts plunging. We thus dub the circular orbit at r = r\ 
the "effective" ISCO. It is important to notice, however, 
that the stable orbits in the "inner" region rs < r < r 2 
do not necessarily have energy and angular momentum 
larger than those of the effective ISCO. For this reason, 
plunging orbits starting at the effective ISCO may hit 
a potential barrier preventing them from reaching the 
compact object, at least in some regions of the parameter 
space (a, q). We will study this situation in detail in the 
next section. 

As q decreases towards q 2l the values of r\ and r 2 ap- 
proach and eventually coincide for q — q 2 . For q < q 2l 
there are no marginally stable orbits in the radial direc- 
tion, i.e. all circular orbits are radially stable. How- 
ever, the marginally stable orbit in the vertical direction, 
r = r3, still exists and marks the position of the ISCO, 
which is therefore determined by the onset of the vertical 
instability. 

In conclusion, for q > q 2 (a) [with q 2 (a) < 0], the ef- 
fective ISCO for quasi-circular inspirals is given by the 
onset of the radial instability, while for q < q 2 (a) it is 
determined by the onset of the vertical instability. Also, 
let us note that the critical value q 2 may be larger than 
—a 2 , that is, the ISCO can be determined by the vertical 
instability even for oblate objects (although since q 2 < 
these objects need to be less oblate than a Kerr BH). 
This situation is summarized in Fig. [T] where we show 
the regions in the (a, q) plane where compact objects are 
prolate/oblate and more prolate/more oblate than a Kerr 
BH, and the regions where the effective ISCO is deter- 
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mined by the vertical/radial instability. 

In the rest of this paper we will present our results for 
MN spacetimes in terms of the standard Boyer-Lindquist 
coordinates (r, #), which are related to the prolate 
spheroidal coordinates (x, y) and the quasi-cylindrical co- 
ordinates (p, z) used in this section by 

p = Vr 2 -2Mr + a 2 M 2 sin (9, 

z = (r-M) cos 0. (2.20) 

and 

r = kx + M , 

cos0 = y. (2.21) 

Also, hereafter we will use coordinates in which the mass 
M of the MN metric is 1. 



III. THE PLUNGE AND THE FORMATION OF 
THICK DISKS INSIDE THE ISCO 

At the inner edge of the thin accretion disk, which 
corresponds to the effective ISCO defined in the previous 
section, the gas is expected to shed and plunge towards 
the central object. To understand the geometry of this 
plunge, it is convenient to plot the region accessible to 
the shedding gas, i.e. the region accessible to particles 
having the ISCO specific energy and angular momentum 
^isco and Lisco- This "plunge" region is defined by 

Kff(£isco,£isco,r,0) >0, (3.1) 

and automatically includes the ISCO (r = risco, & — 
7r/2), but not necessarily a neighborhood of it. More 
specifically, one can have several scenarios: 

1. If the effective ISCO is determined by the onset of 
the radial instability, there are two possibilities: 

(a) The gas plunges directly into the compact ob- 
ject. As in the Kerr case, the "plunge" region 
and therefore the gas may expand above and 
below the equatorial plane forming a sort of 
"plume", but they remain confined near the 
equatorial plane. 

(b) The gas starts to plunge from the ISCO, but 
gets trapped before reaching the compact ob- 
ject. 

2. If the effective ISCO is determined by the onset of 
the vertical instability, there are three possibilities: 

(a) The gas plunges directly into the compact ob- 
ject. However, because the ISCO is stable in 
the radial direction, the gas does not shed on 
the equatorial plane as in the Kerr case, but 
sheds above and below the equatorial plane, 
forming two separate "plumes" . 



(b) The "plunge" region does not contain 
a neighborhood of the ISCO, because 
V e s (Eisco , ^isco , r, 6) is strictly negative 
around it. However, a small perturbation (e.g. 
a small initial velocity) is enough to allow the 
gas to escape the potential well surrounding 
the ISCO, thus shedding above and below the 
equatorial plane and plunging into the com- 
pact object in two separate plumes. Note that 
this case is therefore very similar to scenario 
(2a). 

(c) The gas sheds above and below the equatorial 
plane from the ISCO, but gets trapped before 
reaching the compact object. 

Scenarios (lb) and (2c), where the gas sheds from the 
inner edge of the thin disk but gets trapped before reach- 
ing the compact object, are clearly non- stationary con- 
figurations, because the gas keeps accumulating between 
the object and the ISCO. For the system to settle onto 
a stationary configuration, the gas would have to over- 
flow the potential barrier separating it from the object, 
but this can be achieved only if the gas forms, inside the 
ISCO, a coherent structure that is capable of further dis- 
sipating energy and angular momentum (e.g. through 
viscosity or magnetic fields). 

In this paper we study these "coherent structures" in- 
side the ISCO using the theory of non self-gravitating, 
stationary and axisymmetric thick disks developed in 
Refs. [2U [22] (see also Refs. [5| [23] [24[ ) . We will briefly 
review this formalism in Appendix [Aj but for the present 
discussion it is sufficient to mention that a thick disk is 
completely determined by specifying its inner or outer 
edge and the angular momentum distribution on the 
equatorial plane. More precisely, one needs to assume 
a certain equatorial distribution for the angular momen- 
tum per unit energy t = L/E, and this distribution com- 
pletely determines the distribution £(r, 0) of the angular 
momentum per unit energy outside the equatorial plane 
(see Appendix |A| for details). Assuming in particular a 
power law, we can write 

i{r) = £(r, = tt/2) = ^isco ( ^ , (3.2) 

Vnsco/ 

where iisco = ^isco/^isco an d /3 is a free parameter, 
which needs to be positive in order for the thick disk to 
be stable [23] . 

In our case, we impose that the outer edge of the thick 
disk coincides with the ISCO (so that the thick disk is fed 
by the gas shedding from the thin disk), and we choose 
the parameter j3 such that the thick disk also presents an 
inner shedding point. This requirement completely fixes 
the value of f) to a certain critical value /3 C rit- The gas 
shedding from the inner edge of the thick disk will then 
plunge into the compact object, and the "plunge" region 
accessible to this gas will be described by 

^eff ( -dinner ? -^innen 0)>0, (3.3) 
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1.5 2.0 2.5 3.0 1 2 3 4 5 6 
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FIG. 2. Accretion in Manko-Novikov spacetimes with a = 0.1 and q — 0.2, —0.5, —0.6 and q — —4. The solid black line 
at small radii is the partial horizon x = 1, the region where closed timelike curves exist (i.e. the region where < 0) is 
shown in green/light gray, the red dashed line is the boundary of the ergoregion (i.e. g t t = on that line), while the red 
dot on the equatorial plane marks the position of the ISCO, i.e. the inner edge of the thin accretion disk. The thick disk, if 
it forms, sheds from the ISCO and is denoted by concentric contours, whose label is an upper limit to the log 10 of the gas 
temperature in K (see text for details). In blue/dark grey is the region accessible to the gas shedding from the inner edges of 
the thick disk [i.e. the region where V e s (dinner, dinner, r, 0) > 0, dinner and Lmner being the energy and angular momentum of 
the gas at the inner edges of the thick disk]. If no thick disk is present, in blue/dark grey is the region accessible to the gas 
shedding from the inner edge of the thin disk [i.e. the region where V e s (^isco, £isco, r, 0) > 0]: if this "plunge" region does 
not reach the ISCO, we also show (with solid black lines around the ISCO) the contours V e s (^isco, ^isco, r, 0) = -10- 4 and 
V e ft (Eiscoj ^isco, t, 0) — — 10 -3 , which represent the region where the gas reaching the ISCO can shed if subject to a small 
perturbation. 

For q = 0.2 the ISCO is radially unstable, and the gas plunges directly into the compact object remaining roughly on the 
equatorial plane, as in the Kerr case [scenario (la)]. For q = —0.5 the ISCO is radially unstable and the gas plunges, but does 
not reach the compact object; instead, it gets trapped between the object and the ISCO and forms a thick disk [scenario (lb)]. 
For q = —0.6 the ISCO is vertically unstable, and the gas plunges directly into the compact object outside the equatorial plane 
[scenario (2a)]. For q = — 4 the ISCO is vertically unstable, and the gas does not plunge directly but remains trapped in the 
vicinities of the ISCO [because V e s (£isco, £isco, r, 0) < near the ISCO]. However, the potential barrier is not tall enough to 
allow a thick disk to form, and a small perturbation is enough to cause the gas to fall into the central object [scenario (2b)]. 



where the specific energy and angular momentum of the 
gas at the inner edge of the thick disk, dinner and Li nner , 
are determined by solving simultaneously L- inneT / E- inneT = 

dinner? dinner) and T^gff (-fanner i -^inner i dinner ? dinner) — 0- 

Examples of the scenario outlined above are shown in 
Figs. |2] (for a = 0.1), [3] (for a = 0.5) and |] (for a = 0.9), 
for various values of the anomalous quadrupole moment 



q. The solid black line at small radii is the partial hori- 
zon x = 1, the region where closed timelike curves exist 
(i.e. the region where < 0) is shown in green/light 
gray, while the red dashed line is the boundary of the 
ergoregion (i.e. g u = on that line). The position 
of the ISCO is marked by a red dot on the equatorial 
plane. The thick disk, if it forms, sheds from the ISCO 
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1.5 2.0 2.5 3.0 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 

r sin(0) r sin(0) 

FIG. 3. Same notation as in Fig. [2] but for a = 0.5 and q = 0.2, —0.2, —0.5 and —2. The accretion properties too are 
qualitatively similar, case by case, to those displayed for a = 0.1 in Fig. [2] 



and is denoted by concentric contours, whose label is 
log 10 T, where T is the gas temperature (in K) assum- 
ing a polytropic equation of state p = np^ (po being 
the rest mass density) with T = 2. As we will show 
in Appendix [A[ the temperature for other values of the 
polytropic index T can be obtained by the simple rescal- 
ing T(r) = T(r = 2) x 2(r - l)/r. Because the factor 
2(r — l)/r varies between and 1 when Y varies in its 
allowed range 1 < T < 2, the temperature plotted in our 
figures has to be interpreted as an upper value to the real 
temperature. In blue/dark grey is the "plunge" region ac- 
cessible to the gas shedding from the inner edges of the 
thick disk [i.e. the region where V e R {E- innevi L- innevi r, 0) > 
0]. If no thick disk is present, in blue/dark grey is 
the "plunge" region accessible to the gas shedding from 
the inner edge of the thin disk [i.e. the region where 
Kff(^isco,^isco,^^) > 0]: if this "plunge" region does 
not contain a neighborhood of the IS CO [cf. scenario 
(2b) above], we also show (with solid black lines around 
the ISCO) the contours V G fi (^isco, ^isco, r, 0) = -10 -4 
and Ves(Eisco, Li$co, r,0, ) = -10~ 3 , which represent 
the region where the gas reaching the ISCO can move if 
subject to a small perturbation (e.g. if imparted a small 



initial velocity). 

As can be seen, scenarios (la), (lb) and (2a) happen 
for all the three values of the spin parameter that we 
consider. In particular, although a more detailed scan 
of the parameter plane (a, q) would be needed to draw 
ultimate conclusions, it seems that scenario (la) gener- 
ically takes place when q > <?i(a), scenario (lb) when 
#2 (a) < Q < qi(a): and scenario (2a) when q < q2(a). 
However, for q < q2(a) one can also have scenario (2b) 
- which we see for a = 0.1, q — —4 and a — 0.5, 
q = — 2 - or scenario (2c) - which we see for a — 0.9 and 
q = —0.12. We have not been able to identify a scenario 
(2c) for a = 0.1 or a = 0.5. This is because for these 
values of the spin, the potential barrier preventing the 
gas shedding at the ISCO from plunging above and be- 
low the equatorial plane is smaller than for a = 0.9, and 
a small perturbation is sufficient for the gas to plunge, 
thus making the formation of a stationary thick disk im- 
possible or at least very difficult. Mathematically, thi s 
means that when varying the exponent /3 in Eq. (3.2), 
we were unable to find a value f3 CT \ t for which the thick 
disk has an inner shedding point (and therefore a sta- 
tionary configuration). Scenario (2c) is instead possible 
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FIG. 4. Same notation as in Fig. [2] but for a = 0.9 and q = 0.2, —0.015, —0.02 and —0.12. The accretion properties are 
qualitatively similar to those displayed for a = 0.1 in Fig. [2] and for a — 0.5 in Fig. [3] except for q = —0.12. In that case, the 
ISCO is vertically unstable, and the gas does not plunge directly but remains trapped in the vicinities of the ISCO (because 
V e E (i^isco, £isco, t, 0) < near the ISCO). However, the potential barrier is tall enough to allow the formation of two thick 
disks, above and below the equatorial plane [scenario (2c)]. This contrasts with the cases a = 0.1, q = -4 (Fig. [2} and a = 0.2, 
q — —2 (Fig. [3j), where the potential barrier is smaller and can be easily overcome if the gas is slightly perturbed. 



for a = 0.9 because in that case the potential barrier is 
higher [cf. the contours V<& (^isco, ^isco, r, 0) = -10" 4 
and Kff(£isco,£isco,r,#,) = -10" 3 for a = 0.9 and 
q — —0.12 with those for a = 0.1, q = —4 or a = 0.5, 
q = -2}. 

Finally, let us comment on where the "surface" of the 
compact object may be located. As already mentioned, 
the MN metric is an exact vacuum solution of the Ein- 
stein equations, and describes the spacetime outside a 
generic compact object made of exotic matter, which de- 
termines the free parameters of the MN metric (e.g. the 
anomalous quadrupole moment q). This object should 
of course cover the pathological features of the MN so- 
lutions, such as the closed timelike curves, the curvature 
singularity at x = 1, y = and the partial horizon x = 1, 
\y\ < 1. Even if this were the case, however, this scenario 
would still present some difficulties. First, such an exotic 
compact object would probably be subject to gravita- 
tional instabilities on a dynamical timescale, unless its 
radius is large enough to cover the whole ergoregion. In 



fact, the so-called ergoregion instability is known to hap- 
pen on a dynamical timescale for objects whose exterior 
metric resembles that of Kerr BH, at least for spin pa- 
rameters a < 2 [25] . Second, if the compact object has a 
surface, the kinetic energy of the plunging material must 
eventually be emitted by the surface, while in the case 
of a BH it simply gets lost in the event horizon. This 
argument, coupled with the observation that BH candi- 
dates are indeed dimmer than systems known to contain 
neutron stars, or at least stars with a surface, is a strong 
indication of the existence of event horizons [26]. How- 
ever, we should stress that one can in principle imagine 
models escaping this argument (e.g. a gravastar, where 
the radiation emitted by the surface is so redshifted that 
it becomes observationally negligible [27]). 

Another possibility is that GR breaks down near 
the naked curvature singularity and the closed timelike 
curves of the MN geometry. One can therefore imag- 
ine a completely regular spacetime that is described by 
the MN metric away from the singularity and the closed 



9 



timelike curves, and which does not present these patho- 
logical features any more (see Refs. [20] for some explicit 
ideas in this direction). An advantage of this scenario is 
that the horizon of this "regular" MN spacetime would 
allow these objects to have a dimmer luminosity, in agree- 
ment with the observations [26] , and it may even quench 
the ergoregion instability or make it happen on a non- 
dynamical timescale. 

Of course, the MN spacetimes that we consider, irre- 
spective of their actual physical significance, are also an 
ideal tool to set-up null experiments to test the Kerr BH 
hypothesis, since they reduce exactly to Kerr BHs when 
q = 0. Indeed, any experiment pointing at a value of q 
significantly different from would imply that the object 
under consideration is either a compact object different 
from a BH within GR, or a BH or a compact object in 
a gravity theory different from GR. In this sense, MN 
spacetimes can be used not only to test the BH hypoth- 
esis but also to test the gravity theory itself [28 . 



IV. DISCUSSION 

The scenario (la) discussed in the previous section is, 
as we have mentioned, very similar to the accretion pro- 
cess in a Kerr spacetime. At the ISCO, which is ra- 
dially unstable, the gas plunges into the compact ob- 
ject, remaining roughly on the equatorial plane. Because 
this plunge takes place in a dynamical time and with no 
significant dissipation, the infalling gas is not expected 
to emit a significant amount of radiation, and this is 
a crucial ingredient of the Novikov-Thorne model. Un- 
der this assumption, the accretion luminosity is simply 
L aC c = r]Mc 2 , where r] = 1 — £?isco is the radiative ef- 
ficiency. The evolution of the spin parameter is then 
regulated by [29] 



da 



1 £isco 

d\i\M ~ M Eisco 



2a, 



(4.1) 



neglecting the small correction coming from the radiation 
emitted by the disk and captured by the compact object. 

In the case of a Kerr background, the Novikov- 
Thorne model seems to be confirmed by recent three- 
dimensional general relativistic magnetohydrodynamic 
simulations [16] (see however Ref. [17] for a different con- 
clusion), and can therefore be used to interpret the X-ray 
spectra of BH candidates and estimate their spin [30] . It 
seems therefore plausible that the Novikov-Thorne model 
should work also in MN spacetimes, and in principle one 
can generalize the technique used to estimate the spin pa- 
rameter of BH candidates to measure possible deviations 
from the Kerr geometry (e. g. th e anomalous quadrupole 
moment q) [10]. From Eq. ( |4.1[ ), it also follows that the 
accretion process onto a compact object more oblate than 
a Kerr BH can potentially spin the body up to a > 1 [31] . 

Scenarios (2a) and (2b) are quite similar to scenario 
(la). In these cases the ISCO is vertically unstable, so 
the gas plunges above and below the equatorial plane, 



but this infall still happens on a dynamical time and 
with a negligible emission of radiation. This presumably 
means that the Novikov-Thorne model is a reasonable 
description of the thin disk and t hat the spin of the cen- 
tral object evolves according to (4.1). Also, because the 



shape of the plunge region draws the gas towards the 
rotation axis, outflows and jets might possibly form in 
these scenarios if magnetic fields are present. 

In scenarios (lb) and (2c), the assumption of no ra- 
diation emitted inside the ISCO is not valid, because of 
the presence of the thick disks that we discussed in the 
previous section. As a result, the radiative efficiency is 
not rj = 1 — i^isco like in the Novikov-Thorne model, 
but rj = 1 — dinner, while the spi n evo lution of the central 
object is still described by Eq. ( |4.1| ), but with the ratio 
Lisco/Eisco replaced by L inn JjE inner . The values of 
these quantities are shown in Table [I] for the cases shown 
in Figs. [2] [3] and [4] and as can be seen the corrections 
due to the presence of the thick disks inside the ISCO 
are about 1 — 3% for the ratio L- inneT / E- inneTl and about 
10 — 25% when it comes to rj. 

Moreover, the spectrum of the thick disks inside the 
ISCO is quite different from that of the thin disk. The 
temperature of a thin disk scales like T ~ M -0,25 , M 
being the disk's mass, while the temperature of the thick 
disks inside the ISCO is independent of M (see Ap- 
pendix [A] for details) , like in the case of Bondi accretion 
flows. As can be seen from Figs. [2j [3] and [4j we have 
T - 10 9 - 10 10 K « 0.1 - 1 MeV in the thick disks, 
and at such high temperatures the most efficient cooling 
mechanism is thermal bremsstrahlung, whose emission 
rate is briefly reviewed in Appendix |B) 

In the case of a 10 M compact object, for reasonable 
values of the gas density of the thick disks, e.g. 10 12 
particles/cm 3 , the disk is definitely optically thin. The 
radiation emitted by the thick disk thus scales as M 3 , but 
the flux observable on the Earth is completely negligible: 
considering an object at a distance of 1 kpc, the inten- 
sity of the 7-ray spectrum around 0.1 — 1 MeV is lower 
than 10 -15 7 cm -2 s _1 . For a 10 9 M compact object, 
the thick disk inside the ISCO is not necessarily optically 
thin, but if it is, the 7-ray spectrum of the object may 
include a bump around 0.1 — 1 MeV, with the character- 
istic shape of thermal bremsstrahlung (spectrum almost 
constant till energies ~ T and then exponentially sup- 
pressed). If the object is at 10 Mpc from us, the flux 
on the Earth could be around 1 7 cm -2 s _1 , which is 
potentially observable. For instance, the so-called MeV- 
blazars show this feature [32 . Unfortunately, so far the 
spectrum of active galactic nuclei (AGN) is too poorly 
understood to say anything conclusive, but the presence 
of a thick disk inside the ISCO of the supermassive ob- 
jects in galactic nuclei will likely be testable in the future. 

Finally, let us comment that although the presence of 
thick disks inside the ISCO may be used, at least in prin- 
ciple, as an observational signature of the existence of 
non-Kerr compact objects, the scenarios (lb) and (2c) 
where these disks form seem to arise only in limited re- 
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a q Pcrit 1 — -ElSCO 1 — Sinner ilSCo/SlSCO -dinner/ -dinner 



0.1 -0.5 0.0187 0.069 0.075 3.454 3.426 

0.5 -0.2 0.074 0.093 0.117 3.057 2.990 

0.9 -0.015 0.048 0.167 0.195 2.448 2.418 

0.9 -0.12 0.204 0.310 0.332 1.970 1.921 



TABLE I. The effect of the thick disk on the accretion efficiency and on the spin evolution of the central object, for the cases 
considered in Figs. [2] [3] and ^] that present a thick disk inside the ISCO. 



gions of the (a, q) plane. For instance, as we have men- 
tioned in the previous section, scenario (lb) seems to be 
possible only in a narrow region at the transition between 
scenarios (la) and (2a), while scenario (2c) seems only to 
be possible for sufficiently high spins. 

V. CONCLUSIONS 

The 5 — 2OM compact objects in X-ray binary sys- 
tems and the supermassive bodies in galactic nuclei are 
currently thought to be the Kerr BHs predicted by GR. 
The study of the electromagnetic radiation emitted in the 
accretion process onto these objects can be used to in- 
vestigate their actual nature and therefore test the Kerr 
BH paradigm [TQl [TTJ [13] . In this paper we have stud- 
ied the final stages of accretion, when the gas reaches 
the inner edge of the thin accretion disk, located at the 
ISCO, and plunges into the compact object. We find that 
for non-Kerr compact objects this process is much more 
complicated than in the case of Kerr BH. More specifi- 
cally, depending on the spin a and anomalous quadrupole 
moment q of the compact object we find essentially four 
possible scenarios: 

1. The ISCO is radially unstable, and the gas plunges 
into the compact object remaining roughly on the 
equatorial plane and without emitting significant 
radiation. This is the same scenario as in the Kerr 
case. 

2. The ISCO is radially unstable and the gas plunges, 
but does not reach the compact object. Instead, 
it gets trapped between the object and the ISCO, 
forming a thick disk with T < 10 10 K and emitting 
by thermal bremsstrahlung. 

3. The ISCO is vertically unstable, and the gas 
plunges into the compact object outside the equa- 
torial plane and without emitting significant radi- 
ation. 

4. The ISCO is vertically unstable and the gas 
plunges, but does not reach the compact object. 
Instead, it gets trapped between the object and 
the ISCO and forms two thick disks, above and 
below the equatorial plane. These thick disks have 
T < 10 10 K and emit by thermal bremsstrahlung. 



While the second and the fourth of these scenarios 
seem to happen only for objects more prolate than Kerr, 
and even in that case only in a limited region of the pa- 
rameter space (a,q), they are nevertheless possible and 
may be testable with future data. Our results therefore 
show that care is needed when using the Novikov-Thorne 
model with objects that are different from Kerr BHs, as 
the assumption of negligible radiation emission inside the 
ISCO may not be correct. In particular, an excess of 
emission due to the presence of a thick disk in the region 
inside the ISCO might bias measurements of the spin of 
BH candidates towards high values. It may be interest- 
ing to test this possibility with the spectra of high-spin 
BH candidates such as GRS 1915+105 [33 , when future 
more accurate measurements of the distance to this ob- 
ject will be available. For the time being, however, de- 
spite the peculiar features of the accretion process that 
we discovered in this paper, compact objects more pro- 
late than Kerr BHs cannot be ruled out by astrophysical 
observations. 
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Appendix A: The theory of non self-gravitating, 
stationary and axisymmetric fluid configurations 

Here we briefly review the theory of non self- 
gravitating, stationary and axisymmetric fluid configura- 
tions in generic stationary and axisymmetric background 
spacetimes. For more details, we refer the reader to 
Refs. [2U[22] (but see also Refs. [23l [24] ) . 

We start by considering a perfect fluid with 4-velocity 
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u flmd , which is described by the stress-energy tensor 



PohuZ 



fluid^fluid 



-pg> 



(Al) 



Here p, po? P an d h = (p + p)/po are the pressure, rest- 
mass density, energy density and specific enthalpy. We 
also assume that the fluid is described by a polytropic 
equation of state p = np^ = pqe(T — 1), where £ = p/po~ 
1 is the internal energy per unit rest-mass, and n and V 
are respectively the polytropic constant and polytropic 
index. Because we are neglecting the self-gravity of the 
fluid, g is the background metric (i.e., in our case, the 
MN metric). We assume that the fluid moves on non- 
geodesic circular orbits with 4- velocity 



u fluid = A(r : 0) 



U(r,0)\ 



fl(r, 0) 



d_ 
dt 

dt + £(r, ff)d<j>] 



d 



(A2) 



x fluid 



^ Uld is the energy per 

.fluid /..fluid 



where = u^ uid /u^ uid is the angular velocity, A 
is called the redshift factor, U = 

unit mass as measured at infinity and £ = —u^ md /u^ x 
is the angular momentum per unit energy as measured at 
infinity. Note that £ is conserved for a stationary and ax- 
isymmetric flow in a stationary and axisymmetric space- 
time [34]. The angular momentum per unit energy and 
the angular velocity are related by 



9t<f> + 9ttk 



6« 



9u + gt<t>to 



while u fluid • u fluid = -1 implies 



U 



A 



/ # u £ 2 + 2# t0 £ + 9^ ' 

/ -■ 

9tt + 2#^n + #<^ 2 
AU= 1 



TW 



(A3) 



(A4) 



(A5) 



(A6) 



To calculate the structure of the fluid configuration, 
we use Euler's equation 



^fluid 



(A7) 



where a^ uid is the fluid's 4-acceleration. If the pressure 
is assumed to depend only on r and (which follows 
from the stationarity and axisymmetry of the fluid and 
geometry) and if the equation of sta te is barotropic [i.e., 
if p = p(p)Q one can express Eq. ( |A7[ ) in terms of the 
gradient of a scalar potential W(p): 



..fluid 



W(p) 



I P' 



dp' 



P(p') ' 



(A8) 



Also, from the definition of 4-acceleration and using Eqs. 
flA2] ), ( |A5| , dA6| and the Killing equations V (/i ^) = 
for £ = d/dt and £ = d/d(j), one obtains 



^fluid 



= dfl W = = d„lnU - 



P + P - i-nZ^- (A9) 



Taking the derivative of this equation, anti- 
symmetrizing and using the trivial fact that d^W = 
d^u]^ = d[n U ]U = 0, we obtain d^Qd^i = 0. This im- 
plies V$7 oc V£, and therefore i and Cl have the same 
contour levels [25 • One can th en e xpress $7 as a function 
of £, n = and write Eq. (A9) in the integral form 



W - Wo 



lnC/-lnC/ 



dp' 



o P' + 



P(P') 
n(£')di' 

< out i-JW 



(A10) 



where Wout and £ out are the potential and angular mo- 
mentum per unit energy at the outer edge of the fluid 
configuration. (Of course, W out and £ out can be replaced 
by the values of W and £ at any other point of the fluid's 
boundary, e.g. the inner edge of the configuration if that 
is present.) 

In practice, in order to calculate W — Wout, one spec- 
ifies the functional form of £ on the equatorial pl ane , 
£{r) = £(r, = 7r/2). From that, one can use Eq. (A3) 
to obtain the angular velocity on the equatorial plane, 
Q(r) = (r, = 7r/2), and calculate the integral in 
Eq. ( |A10D as 



1 - {!(£')£' 



ft(r') d£(r') 
I - Ct(r')£(r') dr' 



dr f . (All) 



To calculate In U (r, 0) , one needs instead to reconstruct 
the functional form of £(r, 0) outside the equatorial plane. 
This can be done by solving fi(r, 0) = £l(r) for the equa- 
torial radius f where the angular velocity has the same 
value as at a giv en p oint (r, 0) outside the equatorial 
plane. Using Eq. (A3) and the fact that £l(r, 0) = Q(f) 
implies £{r, 0) — 1(f) (because £ and Q have the same 
contour levels), this equation becomes 



gtgfc 0)+gu(r, 0)£(f) = gt^f) + gttf)£(f) 
g^ (r, 0) + g t <j> (r, 0)£(f) g^ (f ) + g tcf) (f)£(f) 



(A12) 



where we have denoted with g^ v the metric on the equa- 
torial plane. This equation is known as the equation of 
the von Zeipel cylinders, and by allowing one to express 
f as a function of r and it permits calculating both 
$1(r, 0) and £(r,0) outside the equatorial plane. 

The integral Euler equation (A10) can be further sim- 



plified if the equation of state is polytropic, because in 
this case 



1 This is the case for a polytropic equation of state, as p = p/(T - 
1) + (p/,s)i/r 



f p dp' h 

JO P' + Pi?') h out 



(A13) 
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where h out is the specific enthalpy at the outer edge of 
the fluid configuration. Because for a polytropic equation 
of state the enthalpy is 



temperature is independent of the mass of the fluid con- 
figuration, and that the factor 2(T — 1)/T varies between 
and 1 when T varies in its allowed range 1 < T < 2. 



h = l 



1 



Wo 1 



(A14) 



one has h out = 1 (because p = po = Oat the o uter ed ge of 
the fluid configuration), and Eqs. (A10) and (A13) give 



Po(r,0) 



r - 1 r e Wout-w(r,9) _ : j 



i/(r-i) 



, (A15) 



Appendix B: Thermal bremsstrahlung 

The emission rate per unit volume for a fully ionized 
hydrogen gas due to bremsstrahlung processes is [36] 



A, 



brems 



ar 2 m e c s n 2 F(Q e ) , 



(Bl) 



which makes it clear that the fluid occupies the region 
where W(r,0) < W out . 

Once the rest-mass distribution is known, the total rest 



mass is 



M t , 



PoV^4 uid d 3 a; j 



(A16) 



where d 3 x = dr dO d(j) is the coordinate 3- volume element, 
while the mass-energy is 



J (T r r + Ij + Tg - Tt)y/=jdPx . (A17) 



From Eq. ( |A15| | one can also calculate the temperature 
using p — nk B T (where k B is the Boltzmann constant and 
n is the number density). Using the polytropic equation 
of state, one gets 



_ 2(r - 1) 



w c 



>ut -w(r,e) _ -n 



k B r 2k B 

(A18) 

where we have assumed that the mass of the gas parti- 
cles is given by the proton mass ra p . We notice that the 



where a = 1/137 is the electromagnetic fine structure 
constant, r e is the classical electron radius, m e is the elec- 
tron mass, c is the speed of light, n is the number density 
of electrons (protons), and O e = kBT e /rri2C 2 is the di- 
mensionless electron temperature. Here F is the dimen- 
sionless radiation rate due to relativistic bremsstrahlung. 
It is convenient to write F = F ee + F ep , where F ee is the 
contribution from electron-electron collisions, while F ep 
the one from electron-proton collisions. For e < 1, the 
approximated formulas are 



F P , 



20 



(44 - 37T 2 )9^ 5 (1 + l.ie e + O'i - 1.256^) 



2.5\ 



Fe P = 3 ^J-e° e Hi + use*- 34 ) 



For 6 P > 1, we have instead 



(B2) 



F ee = 246 e [ln(2^6 e ) + 1.25] , 

F ep = 12B e [ln(2^6 e + 0.42) + 1.5] , (B3) 

where tje = exp(— 7^) and 7^ ~ 0.5772 is the Euler's 
number. 
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